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Abstract  - For  imaging  and  remote  sensing 
applications,  the  phase  noise  of  the  THz 
generator  must  be  minimized.  This  can  be  done 
as  part  of  the  optimization  process  on  the  basis  of 
a properly  chosen  figure  of  merit.  For  this,  the 
present  paper  provides  an  analytical  phase  noise 
expression  as  a function  of  the  device  parameters 
and  operation  point.  The  quantum  1/f  theory  is 
used  to  calculate  from  first  principles  the  1/f 
noise  present  in  the  device  parameters  and  in  the 
resulting  system  frequency  from  resonant 
tunneling  diodes  (RTDs),  super-electronic  lattice 
devices  (SLEDs),  Gunn  devices  (TEDs),  and 
transit  time  diodes.  In  general,  quantum  /f 
fluctuations  in  the  dissipative  elements  lead  to  a 
Q"4  dependence  of  the  spectral  density  of 
fractional  frequency  fluctuations  and  of  the 
corresponding  phase  noise,  where  Q is  the 
quality  factor. 

I.  INTRODUCTION 

Fluctuations  with  a spectral  density 
proportional  to  1/f  are  found  in  a large  number  of 
systems  in  science,  technology  and  everyday  life. 
These  fluctuations  are  known  as  1/f  noise  in 
general.  They  have  first  been  noticed  by  Johnson  in 
early  amplifiers,  have  limited  the  performance  of 
vacuum  tubes  in  the  thirties  and  forties,  and  have 
later  hampered  the  introduction  of  semiconductor 
devices. 

The  present  paper  is  focused  on  the  general 
origin  of  fundamental  1/f  noise  as  a universal  form 
of  chaos,  and  on  the  cause  of  its  ubiquity.  It  starts 
with  a special  case  of  the  general  1/f  noise 
phenomenon,  the  Quantum  1/f  Effect  (with  its 
conventional  and  coherent  contributions)  which  is 
as  fundamental  as  time  and  space.  The  1/f 
fluctuations  are  a necessary  consequence  of  the 
mathematical  homogeneity  of  the  dynamical  (or 
physical)  equations  describing  the  motion  of  an 
arbitrary  chaotic  or  stochastic  nonlinear  system.  A 
sufficient  criterion  was  derived  by  the  author.  It 
indicates  if  an  arbitrary  system  governed  by  a given 
sys-tem  of  nonlinear  integro-differential  equations 
will  exhibit  1/f  noise.  The  criterion  was  applied  to 
several  particular  systems,  and  is  used  to  predict  the 
fundamental  quantum  1/f  effect  as  a special  case. 


II.  CONVENTIONAL  QUANTUM  1/f  EFFECT 

This  effect1'6  is  present  in  any  cross  section  or 
process  rate  involving  charged  particles  or  cur-rent 
carriers.  The  physical  origin  of  quantum  1/f  noise  is 
easy  to  understand.  Consider  for  example  Coulomb 
scattering  of  current  carriers,  e.g.,  electrons  on  a 
center  of  force.  The  scattered  electrons  reaching  a 
detector  at  a given  angle  away  from  the  direction  of 
the  incident  beam  are  described  by  DeBroglie 
waves  of  a frequency  corresponding  to  their  energy. 
However,  some  of  the  electrons  have  lost  energy  in 
the  scattering  process,  due  to  the  emission  of 
bremsstrahlung.  Therefore,  part  of  the  outgoing 
DeBroglie  waves  is  shifted  to  slightly  lower 
frequencies.  When  we  calculate  the  probability 
density  in  the  scattered  beam,  we  obtain  also  cross 
terms,  linear  both  in  the  part  scattered  with  and 
without  bremsstrahlung.  These  cross  terms  oscillate 
with  the  same  frequency  as  the  frequency  of  the 
emitted  bremsstrahlung  photons.  The  emission  of 
photons  at  all  frequencies  results  therefore  in 
probability  density  fluctuations  at  all  frequencies. 
The  corresponding  current  density  fluctuations  are 
obtained  by  multiplying  the  probability  density 
fluctuations  by  the  velocity  of  the  scattered  current 
carriers.  Finally,  these  current  fluctuations  present 
in  the  scattered  beam  will  be  noticed  at  the  detector 
as  low  frequency  current  fluctuations,  and  will  be 
interpreted  as  fundamental  cross  section  fluctuations 
in  the  scattering  cross  section  of  the  scatterer.  While 
incoming  carriers  may  have  been  Poisson 
distributed,  the  scattered  beam  will  exhibit  super- 
Poissonian  statistics,  or  bunching,  due  to  this  new 
effect  which  we  may  call  quantum  1/f  effect.  The 
quantum  1/f  effect  is  thus  a many-body  or  collective 
effect,  at  least  a two-particle  effect,  best  described 
through  the  two-particle  wave  function  and  two- 
particle  correlation  function. 

Let  us  estimate  the  magnitude  of  the  quantum  1/f 
effect  semiclassically  by  starting  with  the  classical 
(Larmor)  formula  2q2a2/3c3  for  the  power  radiated 
by  a particle  of  charge  q and  acceleration  a.  The 
acceleration  can  be  approximated  by  a delta 
function  a(t)  = Av5(t)  whose  Fourier  transform 
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Av  is  constant  and  is  the  change  in  the  velocity 
vector  of  the  particle  during  the  almost 
instantaneous  scattering  pro-cess.  The  one-sided 
spectral  density  of  the  emitted  bremsstrahlung 
power  4q2(Av)2/3c3  is  therefore  also  constant. 
The  number  4q2(Av)2/3hfc3  of  emitted  photons 
per  unit  frequency  interval  is  obtained  by  dividing 
with  the  energy  hf  of  one  photon.  The  probability 

amplitude  of  photon  emission  [4q(Av)2/3hfc3]1,2e*Y 
is  given  by  the  square  root  of  this  photon  number 

spectrum,  including  also  a phase  factor  e'Y.  Let  i|/ 
be  a representative  Schrodinger  catalogue  wave 
function  of  the  scattered  outgoing  charged  particles, 
which  is  a single-particle  function,  normalized  to 
the  actual  scattered  particle  concentration.  The  beat 
term  in  the  probability  density  p=|vp|2  is  linear  both 
in  this  bremsstrahlung  amplitude  and  in  the  non- 
bremsstrahlung  amplitude.  Its  spectral  density  will 
therefore  be  given  by  the  product  of  the  squared 
probability  amplitude  of  photon  emission 
(proportional  to  1/f)  with  the  squared  non- 
bremsstrahlung  amplitude  which  is  independent  of 
f.  The  resulting  spectral  density  of  fractional 
probability  density  fluctuations  is  obtained  by 
dividing  with  |\|/|4  and  is  therefore 

|v|/|-4SM2(f)  = 8q2(Av)2/3hfNc3  = 

2aA/fN  = r2Sj(f),  (1) 

where  a = e2/h,'  C = 1/137  is  the  fine  structure 
constant  and  aA  = 4q2(Av)2/3hc3  is  known  as  the 
infrared  exponent  in  quantum  field  theory,  and  is 
known  as  the  quantum  1/f  noise  coefficient,  or 
Hooge  constant,  in  electrophysics. 

The  spectral  density  of  current  density  fluctuations 
is  obtained  by  multiplying  the  probability  density 
fluctuation  spectrum  with  the  squared  velocity  of 
the  outgoing  particles.  When  we  calculate  the 
spectral  density  of  fractional  fluctuations  in  the 
scattered  current  j,  the  outgoing  velocity  simplifies, 
and  therefore  Eq.  (1)  also  gives  the  spectrum  of 
current  fluctuations  Sj  (f) , as  indicated  above.  The 

quantum  1/f  noise  contribution  of  each  carrier  is 
independent,  and  therefore  the  quantum  1/f  noise 
from  N carriers  is  N times  larger;  however,  the 
current  j will  also  be  N times  larger,  and  therefore  in 
Eq.  (1)  a factor  N was  included  in  the  denominator 
for  the  case  in  which  the  cross  section  fluctuation  is 
observed  on  N carriers  simultaneously. 

The  fundamental  fluctuations  of  cross  sections  and 
process  rates  are  reflected  in  various  kinetic 
coefficients  in  condensed  matter,  such  as  the 
mobility  p and  the  diffusion  constant  D,  the  surface 
and  bulk  recombination  speeds  s,  and  recombination 
times  x,  the  rate  of  tunneling  jt  and  the  thermal 
diffusivity  in  semiconductors.  Therefore,  the 


spectral  density  of  fractional  fluctuations  in  all  these 
coefficients  is  given  also  by  Eq.  (1). 

When  we  apply  Eq.  (1)  to  a certain  device,  we  first 
need  to  find  out  which  are  the  cross  sections  a or 
process  rates  which  limit  the  current  I through  the 
device,  or  which  determine  any  other  device 
parameter  P,  and  then  we  have  to  determine  both  the 
velocity  change  Av  of  the  scattered  carriers  and  the 
number  N of  carriers  simultaneously  used  to  test 
each  of  these  cross  sections  or  rates.  Then  Eq.  (1) 
provides  the  spectral  density  of  quantum  1/f  cross 
section  or  rate  fluctuations.  These  spectral  densities 
are  multiplied  by  the  squared  partial  derivative 
(dl/da)2  of  the  current,  or  of  the  device  parameter  P 
of  interest,  to  obtain  the  spectral  density  of 
fractional  device  noise  contributions  from  the  cross 
sections  and  rates  considered.  After  doing  this  with 
all  cross  sections  and  process  rates,  we  add  the 
results  and  bring  (factor  out)  the  fine  structure 
constant  a as  a common  factor  in  front.  This  yields 
excellent  agreement  with  the  experiment  in  a large 
variety  of  samples,  devices  and  physical  systems. 

Eq.  (1)  was  derived  in  second  quantization,  using 
the  commutation  rules  for  boson  field  operators.  For 
fermions  one  repeats  the  calculation  replacing  in  the 
derivation  the  commutators  of  field  operators  by 
anticommutators,  which  yields1'6 

p-2Sp(f)  =j-2Sj(f)  =CT_2StT(f)  =2aA/f(N-l)  (2) 

This  causes  no  difficulties,  since  N>2  for  particle 
correlations  to  be  defined,  and  is  practically  the 
same  as  Eq.  (1),  since  usually  N»l.  Eqs.  (1)  and 
(2)  suggest  a new  notion  of  physical  cross  sections 
and  process  rates  which  contain  1/f  noise,  and 
express  a fundamental  law  of  physics,  important  in 
most  high-technology  applications. 

We  turn  now  to  the  connection  to  the  coherent 
Quantum  1/f  Effect,  essentially  caused  by  the 
uncertainty  of  the  electron  mass,  by  the  coherent 
state  of  the  field  of  the  electron.  The  coherent  state 
has  an  uncertain  energy.  The  coherent  state  in  a 
conductor  or  semiconductor  sample  is  the  result  of 
the  experimental  efforts  directed  towards 
establishing  a steady  and  constant  current,  and  is 
therefore  the  state  defined  by  the  collective  motion, 
i.e.  by  the  drift  of  the  current  carriers.  It  is 
expressed  in  the  Hamiltonian  by  the  magnetic 
energy  Em,  per  unit  length,  of  the  current  carried  by 
the  sample.  In  very  small  samples  or  electronic 
devices,  this  magnetic  energy 

Em  = j(B2/87i)d3x  = [nevS/c]2ln(R/r)  (3) 

is  much  smaller  than  the  total  kinetic  energy  E^  of 
the  drift  motion  of  the  individual  carriers 

Ek  = Eimiv2/2  = nSmv2/2  = EJs.  (4) 


Here  we  have  introduced  the  magnetic  field  B,  the 
carrier  concentration  n,  the  cross  sectional  area  S 
and  radius  r of  the  cylindrical  sample  (e.g.,  a current 
carrying  wire),  the  radius  R of  the  electric  circuit, 
and  the  "coherence  ratio" 

s=Em/Ek  =2ne2S/mc2ln  (R/r)  = 2e2N'/mc2,  (5) 

where  N'  =nS  is  the  number  of  carriers  per  unit 
length  of  the  sample  and  the  natural  logarithm 
ln(R/r)  has  been  approximated  by  one  in  the  last 
form.  We  expect  the  observed  spectral  density  of 
the  mobility  fluctuations  to  be  given  by  a relation  of 
the  form 

(l/p2)SM(f)=[l/(l+s)][2aA/fN] 

+[s/(l+s)][2a/7tfN]  (6) 

which  can  be  interpreted  as  an  expression  of  the 
effective  Hooge  constant  if  the  number  N of  carriers 
in  the  (homogeneous)  sample  is  brought  to  the 
numerator  of  the  left  hand  side.  In  this  equation 
aA=2a(Av/c)2/37i  is  the  usual  nonrelativistic 
expression  of  the  infrared  expo-nent,  present  in  the 
familiar  form  of  the  conventional  quantum  1/f 
effect1'6.  This  equation  is  limited  to  quantum  1/f 
mobility  (or  diffusion)  fluctuations,  and  does  not 
include  the  quantum  1/f  noise  in  the  surface  and 
bulk  recombination  cross  sections,  in  the  surface 
and  bulk  trapping  centers,  in  tunneling  and  injection 
processes,  in  emission  or  in  transitions  between  two 
solids. 

Note  that  the  coherence  ratio  s introduced  here 
equals  the  unity  for  the  critical  value  N'  = N"  = 
2*1012/cm.,  e.g.  for  a cross  section  S = 2*10'4  cm2  of 
the  sample  when  n = 1016.  For  small  samples  with 
N'«N"  only  the  first  term  survives,  while  for 
N'>N"  the  second  term  in  Eq.  (6)  is  dominant. 

III.  DERIVATION  OF  MOBILITY  QUANTUM  1/f 
NOISE  IN  N+-P  DIODES 

For  a diffusion  limited  n+-p  junction  the  current  is 
controlled  by  diffusion  of  electrons  into  the  p - 
region  over  a distance  of  the  order  of  the  diffusion 
length  L = (Dnxn)1/2  which  is  shorter  than  the  length 
wp  of  the  p - region  in  the  case  of  a long  diode.  If 
N(x)  is  the  number  of  electrons  per  unit  length  and 
Dn  their  diffusion  constant,  the  electron  current  at  x 
is 

Ind  = - eDndN/dx,  (7) 

where  we  have  assumed  a planar  junction  and  taken 
the  origin  x = 0 in  the  junction  plane.  Diffusion 
constant  fluctuations,  given  by  kT/e  times  the 
mobility  fluctuations,  will  lead  to  lo-cal  current 
fluctuations  in  the  interval  Ax 

8AInd(x,t)  = IndAx8Dn(x,t)/Dn. 


The  normalized  weight  with  which  these  local 
fluctuations  representative  of  the  interval  Ax 
contribute  to  the  total  current  Id  through  the  diode  at 
x = 0 is  determined  by  the  appropriate  Green 
function  and  can  be  shown  to  be  (l/L)exp(-x/L) 
for  wp/L  » 1.  Therefore  the  contribution  of  the 
section  Ax  is 

5AId(x,t) 

=(Ax/L)exp(-x/L)Ind8D„(x,t)/Dn,  (9) 

with  the  spectral  density 
SAld(x>f) 

=(Ax/L)2  exp(-2x/L)  Ind2  SDn(x,f)/Dn2  (10) 

For  mobility  and  diffusion  fluctuations  the 
fractional  spectral  density  is  given  by  aHnd/fNAx, 
where  aHnd  is  determined  from  quantum  1/f  theory 
according  to  Sec  II.  With  Eq.  (7)  we  obtain  then 

SAid(x,f) 

= (Ax/L2)  exp  (-2x/L)  (eDndN/dx)2aHnd/fN  (11) 

The  electrons  are  distributed  according  to  the 
solution  of  the  diffusion  equation,  i.e. 

N(x)  = [N(0)  - Np]  exp(-x/L); 
dN/dx  = -{[N(0)-Np]/L}  exp(-x/L).  (12) 

Substituting  into  Eq.  (11)  and  simply  summing  over 
the  uncorrelated  contributions  of  all  intervals  Ax, 
we  obtain 

SId(f)  = a„nd(eDn/L2)2J[N0-Np]2e-4x/Ldx  / 

{[N(0)  -Np]  e‘x/L  + Np}  (13) 

The  integral  is  from  0 to  w.  We  note  that 
eD„/L2  = e/xn.  With  the  expression  of  the  saturation 
current  I0  = e(Dn/xn)1/2Np  and  of  the  current 
I = I0[exp(eV/kT)  - 1],  we  can  carry  out  the 
integration  from  0 to  1 

Sid(f)  = aHnd(eI/fxn)  J du/(au  + 1) 

= aHnd(eI/fxn)  F(a).  (14) 

Here  we  have  introduced  the  notations 

u = exp(-x/L),  a = exp(eV/kT)  - 1, 

F(a)  = l/3-l/2a  + 1/a2  - (1/a3)  ln(l+a)  (15) 

Eq.  (14)  gives  the  diffusion  noise  as  a function  of 
the  quantum  1/f  noise  parameter  aHnd.  A similar 
result  can  be  derived  for  the  quantum  1/f 
fluctuations  of  the  recombination  rate  r in  the  bulk 
of  the  p - region,  the  only  difference  being  the 
presence  of  aHnr  instead  of  aHnd  in  Eq.  (14).  The 
total  noise  is  the  given  by  Eq.  (14)  with  aHnd 
replaced  by  the  sum  aHnd  + aHnr 

SId(f)  = (a„nd  + aHnr)(eI/fxn)  F(a).  (16) 
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IV.  1/f  FLUCTUATIONS  IN  RTDs 

Resonant  tunneling  diodes  have  been  proposed  as 
generators  of  THz  oscillations  and  radiation.  They 
consist  of  two  potential  barriers  enclosing  a 
quantum  well.  Electrons  penetrating  the  potential 
barriers  by  tunneling  are  controlled  by  the 
quasistationary  energy  levels  defined  by  the 
potential  well.  If  their  energy  is  close  to  the  first 
energy  level  in  the  well,  resonance  occurs,  and  a 
peak  Ip  of  the  current  through  the  diode  occurs. 
This  corresponds  to  an  applied  bias  voltage  VP.  If, 
however,  the  applied  voltage  increases  further,  only 
a negligibly  small  non-resonant  current  trickle 
remains  at  the  voltage  V=Vv  and  a broad  valley  is 
observed  in  the  I/V  characteristic.  Scattering 
processes  that  reduce  the  energy  of  the  carriers  to  a 
value  close  to  eVP  will  always  be  present, 
generating  a finite  current  minimum  Iv  at  Vv. 
Between  VP  and  Vv  there  is  a negative  differential 
conductance 

G=-(IP-Iv)/(Vv-VP)  (17) 

on  the  I/V  curve,  that  is  used  to  generate 
oscillations. 

1/f  noise  in  Iv  is  given  by  Eqs.  (l)-(2)  with 

(Av/c)2=2eVv/m.  (18) 

Taking  for  instance  VP=0.4  V,  IP=2.5  108A/m2, 
Vv=0.6  V,  Iy=4  107  A/m2,  we  obtain  with 
meff=0.068  m0, 

NIV'2  SJv(f)  =2aA/f  =7.4/0.068  10'9  =1.3  10‘7  (19) 

N is  given  by 

N =tlv/e,  (20) 

where  x is  the  life  time  of  the  carriers.  With  a cross- 
sectional  area  of  10‘6cm2  and  t=10‘10,  we  obtain 
N=  2.5  106. 

Finally,  the  quantum  1/f  frequency  fluctuations  can 
be  obtained  from  the  formula 

co2=-G3Lc/Cd3  (Rl+Rs)  , (21) 

which  is  obtained  by  equating  the  RF  circuit-limited 
output  power  with  the  DC  bias  circuit-limited 
power.  Here  Lc  is  the  inductance  in  the  bias  circuit, 
Cp>  is  the  diode  capacitance,  Rp  the  resistance 
introduced  by  the  RF  circuit  inductance,  and  Rs  is 
the  diode  series  resistance.  From  Eq.  (17)  the 
fluctuations  in  G caused  by  Iy  are  calculated  and 
substituted  into  the  linearized  form  of  Eq.  (21). 

V.  1/f  FREQUENCY  FLUCTUATIONS  FROM  1/f 
NOISE  IN  RF  DISSIPATIVE  ELEMENTS 


oscillator  with  losses 

dx/dt  + ydx/dt  + co02x  = F(t)  (22) 

The  quantum  1/f  fluctuations  are  present  in  the  loss 
coefficient  y.  They  are  given  by  an  expression  of 
the  form 

S6y/Y(f)  = A/f  (23) 

where  A is  a quantum  1/f  coefficient  characterizing 
the  elementary  loss  process. 

The  resonance  frequency  is  given  by 

cor2  = 0)o2  + f (24) 

The  quantum  1/f  fluctuations  in  the  resonance 
frequency  are  given  by  cor8cor  = -2ySy,  or 

8co/cor  =-2(y/cor)  28y/y  = ~(1/2Q2)  (Sy/y) 

(25) 

where  Q=cor/2y  is  the  quality  factor.  Squaring, 
averaging  and  particularizing  Eq.  (25)  for  the  unit 
frequency  interval,  we  obtain  the  Q'4  law 

Sy(f)  =<(8tor/cor)  2>f  =(l/4Q4)<(Sy/y)  2>f 

=A/4fQ4,  (26) 

where  y is  the  fractional  frequency  fluctuation 
Scor/cor.  This  law  is  also  applicable  to  quartz 
resonators5.  The  examples  in  the  last  two  sections 
show  how  the  phase  noise  is  to  be  found  in  other 
solid  state  sorces  of  THz  radiation. 
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